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Abstract 

For an elliptic problém with two space dimensions, we propose to formulate the finite 
volume method with the help of Petrov-Galerkin mixed finite elementsthat are based on 
the building of a dual Raviart-Thomas basis. 

Résumé 

Pour un probléme elliptique bidimensionnel, nous proposons de formuler la méthode des 
volumes finiš avec des éléments finiš mixtes de Petrov-Galerkin qui s'appuient sur la 
construction ďune base duale de Raviart-Thomas. 

Keywords Finite volumes, mixed finite elements, Petrov-Galerkin variational formula- 
tion, inf-sup condition, Poisson equation. 



1. Introduction 

• Let Í7 be a bidimensional bounded convex domain in with a polygonal boundary 
dVí. We consider the homogeneous Dirichlet problém for the Lapiace operátor in the 
domain Í7 : 

,^ , í —Au = f in íl 

(1) <^ 

\ u =0 on the boundary dfl of Q. 
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We suppose that the dátum / belongs to the space L^(fž). We introduce the momentum 
p defined by 

(2) p = ^u. 

Taking the divergence of both terms arising in equation ([2]), taking into account the 
relation ([T]), we observe that the divergence of momentum p belongs to the space L^(fž). 
For this reason, we introduce the vectorial Sobolev space ií(div, fž) = {q E L^(Í7) x 
L^(fž) , div g G L^(íž) }. The variational formulation of the problém ([T]) with the help 
of the pair ^ = {u, p) is obtained by testing the definition ([2]) against a vector valued 
function q and integrating by parts. With the help of the boundary condition, it comes 
■ (í*) q) + (^1 div g) = O, Vg G ií(div, Q). Independently, the relations ([T]) and ([2]) are 
integrated on the domain íl after multiplying by a scalar valued function v G L^(íž) . 
We obtain : (divjo, v) + (/, = O, Vf G L^(fž). The "mixed" variational formulation 
is obtained by introducing the product space V defined as V = L^(íž) x H{dw, Q), 
II {u, p) \\y = II u IIq + II p IIq + II divp IIq, the foUowing bilinear form 7(», •) defined 
on V x V : 

(3) 7((m, p), {v, g)) = {p, g) + {u, div g) + (divp, v) 

and the linear form cr(») defined on V according to : < a, ( > = —{f, v), , C = ("í^; ^ 
V. Then the Dirichlet problém ([T]) takes the form : 

(4) I ^ ^ 

Due to classical inf-sup conditions introduced by Babuška |Ba71] . the problém admits 
a unique solution ^ E V . 



2. Mixed Finite Elements 

• We introduce a mesh T that is a bidimensional cellular complex composed in our case 
by triangular elements K {K G S-j-), straight edges a (a G At) and ponctual nodes 
S {S E Sf)- We conside also classical finite dimensional spaces L^(fž) and Hf{dw, Q) 
that approximate respectively the spaces L^(fž) and ií(div, Q). A scalar valued function 
v G L^(Í7) is constant in each triangle K of the mesh : L^(fž) = |f : Í7 — > M, Vií G 
S-j-i G M, V a; G K, v{x) = vk} ■ vector valued function function g G H-j-{div, Q) 
is a linear combination of Raviart-Thomas |RT77] basis functions oí lower degree, 
defined for each edge a G Ar as follows. 

• Let a G At be an internal edge of the mesh, denote by S and the two vertices 
that compose its boundary (see Figúre 1) : da = { S, N } and by K and L the two 
elements that compose its co-boundary d'^a = {K, L} in such a way that the normál 
direction n is oriented from K towards L and that the pair of vectors {n, ŠŤ^) is 
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direct, as shown on Figúre 1. We denote by W (respectively by E) the third vertex of 
the triangle K (respectively of the triangle L) : K = {S, N, W), L = {N, S, E). 




Figúre 1. Co-boundary {K, L) of the edge a = {S, N). 

The vector valued Raviart-Thomas |RT77] basis function (fa is defined by the relations 
^aix) = ^^^{x — W)whenx E K, ipa{x) = — {x — E) when x E L and ipa{x) = O 
elsewhere. When the edge a is on the boundary díl, we suppose that the normál n 
points towards the exterior of the domain, so the element L is absent. We have in all 
cases the iíí(div, íl) conformity : (pa ^ H{dw, íl) and the degrees of freedom are the 
fluxes of vector field for all the edges of the mesh (see |RT77j ) : d7 = 

Sa,b, "^a, b E At- a vector valued function q G //7-(div, Vt) is a linear combination of 
the basis functions Lpa '■ 1 = Sag^r '^'^ ^ H-T-{div, Q) = < ípi, , b E Ar > ■ 
• The mixed finite element method consist in choosing as discrete linear space the 
following product : Vj- = L^(íž) x H-i-{dw, Q) and to replace the letter V hy Vj- 
inside the variational formulation : ^ ^) li^Ti C) =< o", C >) E Vf. In 
other terms 

Íur E Ľij-iVí) , pr G iír(div, íž) 
br, q) + (^r, div g) = 0, Vg G /ír(div, n) 
(divpr,t^) + U.V) =0, ELliO). 

The numerical analysis of the relations between the continuous problém and the dis- 
crete problém ([5]) as the mesh T is more and more refined is classical |RT77) . The above 
method is popular in the context of petroleum and nuclear Industries but suffers from the 
fact that the associated linear systém is quite difíicult to solve from a practical point of 
view. The introduction of supplementary Lagrange multipliers by Brezzi, Douglas and 
Marini |BDM85] allows a simplification of these algebraic aspects, and their interpretation 
by Croisille in the context of box schemes |Cr2k] gives a good mathematical foundation 
of a popular numerical method. 



3. Finite Volumes 

• From a theoretical and practical point of view, the resolution of the linear systém ([5]) 
can be conducted as foUows. We introduce the mass-matrix Ma^h = {'-Pa-, '^b)i o, b E Aj- 
associated with the Raviart-Thomas vector valued functions. Then the first equation of 
([5]) determines the momentum pj- = YlaeA PT,a^a as a function of the mean values 
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(6) PT,a = - ^ i^^'^)a,b / diví/?bdx. 

6e^r ' K££r 

The representation ([H]) suffers at our opinion form a major defect : due to the fact that the 
matrix is fuU, the discrete gradient pr is a global function of the mean values uj^k 

and this property contradicts the mathematical foundations of the derivation operátor to 
be linear and local. An a posteriori correction of this defect has been proposed by 
Baranger, Maitre and Oudin |BM096) : with an appropriate numerical integration of the 
mass matrix M, it is possible to lump it and the discrete gradient in the direction n of 
the edge a is represented by a formula of the type : 

(7) PT,a = 7 

with the notations of Figúre 1. The substitution of the relation ([7]) inside the second 
equation of the formulation conducts to a variant of the so-called finite volume method. 
In an analogous manner, the family of finite volume schemes proposed by Herbin |He95] 
suppose a priori that the discrete gradient in the normál direction admits a representation 
of the form ([7]). Nevertheless, the intuition is not correctly satisfied by a schéme such that 
d?]). The finite difference "'^'^^ "'^^ wish to be a a good approximation of the gradient 
Pj- = Vtĺr in the direction kŽ whereas the coefficient p-j^a is an approximation of 
J^Vu-fU dr in the normál direction (see again the Figúre 1). When the mesh T 
is composed by generál triangles, this approximation is not completely satisfactory and 
contains a real limitation of these variants of the finite volume method at our opinion. 
• In fact, the finite volume method for the approximation of the diffusion operátor 
has been first proposed from empirical considerations. Following e. g. Noh |No64j and 
Patankar |Pa80| . the idea is to represent the normál interface gradient J^Vu-fU dr as 
a function of neighbouring values. Given an edge a, a vicinity V(a) is first determined 
in order to represent the normál gradient pj^a = /a ^u-j-»n dr with a formula of the type 



/ Vur-n dr = V" ga,K ut,k 



Then the conservation equation div p + / = O is integrated inside each celí K E Sj- is 
order to determine an equation for the mean values u^^k for all K G S-j-. The difficulties 
of such approches have been presented by Kershaw |Ke81] and a variant of such schéme 
has been first analysed by Coudiére, Vila and Villedieu |CVV99] . The key remark that 



we have done with F. Arnoux (see |Du89) ). also observed by Faille, Gallouét and Herbin 
|FGH91] is that the representation ([8]) must be exact for linear functions uj- . We took 
this remark as a starting point for our tridimensional finite volume schéme proposed in 
|Du92) . It is also an essential hypothesis for the result proposed by Coudiére, Vila and 
Villedieu. 
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4. Finite volumes as mixed Petrov-Galerkin finite elements 

• In this contribution, we propose to discretize the variational problém (jl]) with the 
Petrov-Galerkin mixed finite element method, first introduced by Thomas and Trujillo 
[TT99]. In the way we have proposed in [Du2kj . the idea is to construct a discrete 
functional space H-j- {div, íl) generated by vectorial functions ip*, a & At, that are 
conforming in the space ií(div, íl) : <p*a ^ ií(div, íž) and to represent exactly the 
dual basis of the family 6 G Af} with the scalar product : ) = 

^a,bi a, b E At- Then i/7- (div, íl) = < ipl , b E Aj- > . The mixed Petrov-Galerkin 
mixed finite element method consists just in replacing the space iÍ7-(div, íl) by the dual 
space iir (div. Q) for test functions in the first equation of discrete formulation ([5]). 
We obtain by doing this the so-called Petrov-Galerkin finite volume schéme : 

Íur G L?j-{Q) , pr G iir(div. Q) 
(PT, q) + (ur, divq) = 0, Vg G iir(div, íl) 
{dwpr,v) + {f,v) =0, eL^íl). 

We introduce a compact form of the previous mixed Petrov-Galerkin formulation with the 
help of the product space definedby = L^íl) x iir (div, ^)- Then the discrete 
variational formulation (|9]) admits the form : ^r ^ ^Ti i{Ít^ C) = <^)C>)^C ^ ■ 



5. Stability analysis 

• We suppose in the following that the mesh T is a bidimensional cellular complex 
composed by triangles as proposed in the previous sections. Following the work of Ciarlet 
and Raviart |CR72] . for any element K E Sj- we denote by the diameter of the 
triangle K and by px the diameter of the inscripted ball inside K. We suppose that 
the mesh T belongs to a family U0 {6 > 0) of meshes defined by the condition T G 
Uq <í==^ y k e £t 1 — < 0. We suppose also that the dual space iir (div, ^) 
constructed by the previous conditions satisfies the following hypothesis. 
Hypothesis 1. Interpolation operátor iir (div, ^) — ^ iir (div, ^) • We sup- 
pose that the mesh T belongs to the family Uq. Let iir(div, íl) 3 q 1 — líg G 
iir (div, ^) be the mapping defined by the condition 

(10) qaiPa) = ^ qaVt, ^ qa & iir(div, íl). 

We suppose that the dual basis is constructed in such a way that there exists strictly 
positive constants A, B, D, E that only depends on the parameter 9 such that we have 
the following estimations : 

^ A\\q\\l < iq,Uq) , Vg G iir(div, fi) 

I II ng llo < S II g llo , Vg G iir(div, íl) 

1 II divHg llo < D \\ divq ||o , Vg G iir(div, íl) 

^ (divg, divHg) > E || divg , Vg G iir (div, fi) . 
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Proposition 1. Technical lemma about lifting of scalar fields. Let 6* be a 

strictly positive parameter. We suppose that the dual Raviart-Thomas basis satisfies the 
Hypothesis 1. Then there exists some strictly positive constant F that only depends on 
the parameter 6 such that for any mesh T that belongs to the family Ug, and for any 
scalar field u constant in each element K E £t (m € L^(íž)), there exists some vector 
field q G ifrídiv, íž) such that || q \\H{án,Q.) ^ F \\ u ||o and {u, div g) > || m ||o . 
Proposition 2 Discrete stability. Let ^ be a strictly positive parameter. We 
suppose that the dual Raviart-Thomas basis satisfies the Hypothesis 1. Then we have 
the foUowing discrete stability for the Petrov- Galerkin mixed formulation ^ : 3 /3 > 
O, vr G W,, Ve e Vr such that || i \\y = 1, G F^, || C \\v ^ 1, 7(^,0 > Z^, 
with 7(», •) defined at the relation ([3]) and /3 chosen such that ^1 — — > 

Theorem 1 Optimal error estimate. Let íž be a two-dimensional open convex 
domain of with a polygonal boundary, u G ií^(íž) be the solution of the problém 
([T]) considered under variational formulation and p = Vu be the associated momentum. 
Let ^ be a strictly positive parameter and Ue a family of meshes T that satisfy the 
Hypothesis L Let ^ = {u-j-, Pt) G V7- be the solution of the discrete problém ([9]). Then 
there exists some constant C > O that only depends on the parameter 9 such that 
\\u-ur llo + \\p-PT \\H{div,n) < C hr \\ f ||o • 



6. Towards a first Petrov- Galerkin finite volume schéme 




Figúre 2. Support V(SN) of the dual Raviart-Thomas basis function (p^j^t- 

Theorem 2 We suppose that the internal edge a links the two vertices S and N 
(see the Figúre where a = SN, O is the middle of SN and n is the associated normál 
direction), if the support of the dual Raviart-Thomas basis function í^gj,^ is the vicinity 
V(a) = {K, L, M, P, Q, R} of the edge a composed by six triangles presented on Figúre 
2 and if the divergence of the dual Raviart-Thomas basis function is equal to a constant 
field in each triangle of V(a) (div(y9*f>j G L^(fž)), then the five mean flux values rj = 
/sN V?SN • ^ dr, a = J^^ ip'l^ • uen dr , /3 = J^^ cp^is • "wn dr , 7 = /^g í^f^ • ^ws dr 
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and S 



(12) 



IsE V^SN^'^SEclr satisfy the following three scalar constraints : 

rjKĹ + aĽÚ + (3vá + -fqŔ + ôlA = |Srt| n 



(13) alÚ.wX + (3FŔ.Eň + -fqŔ.EČ + 5Ľň.W^ = -3 | ŠŤ^ | n,{ÔÍ + (yŔ). 



• The finite volume approach is then obtained in the spirit of ([8]) with a six point 
schéme for the mean gradient in the normál direction thanks to the first equation of the 
mixed variational formulation : 



We remark that the constraints 0121) express that the relation (ÍT^ is exact if the field Uj- 
is an affine function. Taking into account the fact that we have five parameters for the 
definition of the finite volume schéme (relation (ÍT4l) ) and only three constraints (relations 
( ÍT2|) and ( 1T3|) ) for these parameters, the stability seems a reasonable goal, even if the 
problém remains essentially open for generál triangular meshes. 
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